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Optical traps and lattices provide a new opportunity to study strongly correlated high 
spin systems with cold atoms. In this article, we review the recent progress on the 
hidden symmetry properties in the simplest high spin fermionic systems with hyperfine 
spin F = 3/2, which may be realized with atoms of 132 Cs, 9 Be, 135 Ba, 137 Ba, and 
201 Hg. A generic 50(5) or isomorphically, 5p(4)) symmetry is proved in such systems 
with the s-wave scattering interactions in optical traps, or with the on-site Hubbard 
interactions in optical lattices. Various important features from this high symmetry are 
studied in the Fermi liquid theory, the mean field phase diagram, and the sign problem 
in quantum Monte-Carlo simulations. In the s-wave quintet Cooper pairing phase, the 
half-quantum vortex exhibits the global analogue of the Alice string and non-Abelian 
Cheshire charge properties in gauge theories. The existence of the quartetting phase, a 
four-fermion counterpart of the Cooper pairing phase, and its competition with other 
orders are studied in one dimensional spin-3/2 systems. We also show that counter- 
intuitively quantum fluctuations in spin-3/2 magnetic systems are even stronger than 
those in spin-1/2 systems. 

Keywords: high spin; hidden symmetry; quintet Cooper pairing; quartetting. 



1. Introduction 

The past decade has witnessed the great progress in the cold atomic physics. The 
first generation of Bose-Einstein condensation (BEC) in alkali atoms was realized in 
mag netic traps. In such systems, spins of atoms are fully polarized by the Zeemann 
field ^21, thus atoms are effectively of single component. A few years later, impor- 
tant achievement was made to release the spin degrees of freedo m by using optical 
methods to trap atoms, such as optical traps and lattices 3 4 o additional ad- 
vantage of optical lattices is the excellent controllability of the interaction strength 
from the weak coupling to strong coupling regimes. For example, the superfluid to 
Mott insulator transition of bosons has been experimentally observed «k All these 
progresses provide a controllable way to study high spin physics, in particular, the 
strongly correlated high spin physics. 

The high spin physics with cold atoms contains novel features which does not 
appear in solid state systems. In many transition metal oxides, several electrons are 
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combined by Hund's rule to form an onsite composite high spin object with S > h. 
However, the intersite coupling is still dominated by the exchange process of a single 
pair of electrons, thus the leading order spin exchange is the bilinear Heisenberg 
term. Due to the 1/5 effect, quantum fluctuations are typically weak. In other 
words, the solid state high spin systems are more classical than spin-i systems. 
In contrast, such restriction does not happen in high spin systems with cold atoms 
because the building block of such systems, each atom, carries a high hyperfine spin 
F which includes both electron and nuclear spins. This feature renders the possible 
high hidden symmetries and strong quantum fluctuations. 

The high spin bosonic systems exhibit much richer phase diagram than the 
spinless boson. In spin-1 systems of 23 Na and 87 Rb atoms, polar and ferromagnetic 
spinor condensations and related collective modes were studied by Ho ■ and Ohmi 
et al®. The interplay between the U(l) charge and SU(2) spin degrees of freedom 
in the spinor condensate is characterized by a Zi gauge symmetry. The topolo gical 
aspect of this hidden Zi gauge symmetry was investigated extensively by Zhou^iI3 
and Demler et a/OH. Furthermore, the spin si ngle t state and spin nematic state were 
also studied in optical lattices by Zhou et al. and Imambekov et a/ ^jl Recently, 
the spin-2 hyperfine state of 87 Rb atom H^USI and spm _ 3 atom of 52 Cr | 16 | 17 | lg fll 
have attracted a lot of attention. Various spinor condensates and spin ordered Mott 
insulating states have been classified. Remarkably, a biaxial spin-nematic state can 
be stabilized in 52 Cr systems which can support interesting topological defects of 
non-Abelian vortices. 

On the other hand, relatively less works have been done for high spin fermions. 
However, such systems indeed have i nter esting properti es, and thus deserve atten- 
tion. Pioneering works by Yip et al. ^ an d Ho et al. showed that high spin 
fcrmionic systems exhibit richer structures of collective modes in the Fermi liq- 
uid theory, and more diversities in Co ope r pairing patterns. More recently, a large 
progress has been made by Wu et al. on the symmetry properties in spin-3/2 
systems, which may be realized with atoms such as 132 Cs, 9 Be, 135 Ba, 137 Ba, 201 Hg. 
Such systems are very special in that in additional to the obvious spin SU (2) sym- 
metry, they possess a hidden and generic 50(5) or isomorphically, Sp{A) symmetry. 
In this paper, we will review this progress. 

Before we move on, let us briefly review some group theory knowledge. The 
SO (5) group describes the rotation in the 5-dimensional space spanned by five real 
axes ni ~ n 5 which form the vector representation of 50(5). SO (5) group has ten 
generators L a {, (1 < a < b < 5), each of which is responsible for the rotation in 
the ab plane. These generators form 5*0(5) 's 10 dimensional adjoint representation. 
50(5) 's fundamental spinor representation is 4-dimensional. Rigorously speaking, 
its spinor representation is faithful for 50(5) 's covering group 5p(4). The rela- 
tion between 50(5) and Sp{4) is similar to that between 50(3) and SU{2). For 
simplicity, we will not distinguish the difference between SO (5) and Sp(i) below. 

The origin of the 50(5) symmetry in spin-3/2 systems can be explained as 
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follows: The kinetic energy term shows an obvious ST7(4) symmetry because of 
the equivalence among the four spin components ^ a (a = ±|,±i). However, this 
SU(4) symmetry generally speaking is broken explicitly down to 50(5) by inter- 
actions. Pauli's exclusion principle requires that in the s-wave scattering channel, 
only interactions go in the total spin singlet channel (St = 0) and gi in the quintet 
channel (St — 2) are allowed, while those in the channels of St = 1,3 are forbidden. 
Remarkably, the singlet and quintet channels automatically form the identity and 
5-vector representations for an 5*0(5) group (see Sect. [5] for detail), thus making 
the system 50(5) invariant without fine tuning. The same reasoning also applies 
to the existence of the 5*0(5) symmetry in the lattice Hubbard model. The validity 
of this 50(5) symmetry is regardless of dimensionality, lattice geometry and im- 
purity potentials. Basically, it plays the same role of the SU(2) symmetry in the 
spin-^ systems. This 50(5) symmetry purely lies in the particle-hole channel as 
an extension of the spin SU(2) algebra. It is qualitatively different from the 50(5) 
theory in the high T c superconductivity which involves both particle-hole and 
particle-particle channels. 

To our knowledge, such a high symmetry without fine tuning is rare in condensed 
matter and cold atomic systems, and thus it is worthwhile for further exploration. 
Below we outline several important consequences from the 50(5) symmetry and 
other interesting properties in spin-3/2 systems. The 50(5) symmetry brings hidden 
degeneracy in the collective modes in the Fermi liquid theory This symmetry 
greatly facilitates the understanding of the mean field phase diagram in the lattice 
system I^H. Furthermore, due to the time-reversal properties of the 50(5) algebra, 
the sign problem of the quantum Monte-Carlo algorithm is proved to be absent in 
a large part of the phase diagram 1 2 1 1 2 3 1_ 

Spin-3/2 systems can support the quintet Cooper pairing phase i.e., Cooper pair 
with total spin Stot = 2. This high spin superfluid state possesses the topological 
defect of the half-quantum vortex often called the Alice string. An interesting prop- 
erty of the half quantum vortex loop or pair is that they can carry spin quantum 
number as a global analogy of the Cheshire charge in the gauge theories. Interest- 
ingly, in the quintet pairing state, when a quasiparticle carrying spin penetrates 
the Cheshire charged half-quantum vortex loop, quantum entanglement is gener- 
ated between them. The Alice string and non-Abelian Cheshire charge behavior in 
spin-3/2 systems were studied by Wu et a^^H. 

Fermionic systems with multiple components can support multiple-particle clus- 
tering instabil ities, which m eans more than two fermions come together to form 
bound states E 5 | 26|27 | 28 | 29 | 30 | 31 ] p or ex ample, baryons are bound states of three- 
quarks, a-particles are bound states of two protons and two neutrons, and bi- 
excitons are bound state of two electrons and two holes. Spin-3/2 systems can ex- 
hibit quartetting order as a four-fermion counterpart of the Cooper pairing. Taking 
into accoun t the arrival of the fermion pairing superfluidity by Feshbach resonances 
| 32 | 33 | 34 135J ^ j g na ^ ura i i expect the quartetting order as a possible research fo- 
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cus in the near future. The existence of the quartetting order and its competition 

IOOI 

with the singlet pairing order were investigated by Wu — in ID spin-3/2 systems. 
The generalization of the quartetting order to even higher spin systems at ID was 



investigated by Lecheminant et a/ 1 ^. 

Spin-3/2 magnetic systems are characterized by strong quantum fluctuations. 
This is a little bit counter-intuitive because one would expect weak quantum fluc- 
tuations due to the high spin. However, because of the high symmetry of 5*0(5), 
quantum fluctuations are actually even stronger than those in spin ^ systems. For 
example, spin-3/2 systems at quarter filling (one particle per site) can exhibit the 
SU(A) plaquette order which is a four-site counterpart of the dimer state in spin-i 
systems. An SU (4) Majumdar-Ghosh model was constructed by Chen et al. 
in the spin-3/2 two-leg ladder systems. Their ground state is exactly solvable and 
exhibits such plaquette order. 

On the experimental side, in addition to alkali atoms, considerable progress has 
been made in trapping and cooling the divalent alkaline-earth atoms 12 B | 39 | 40 | 4l ] 
Two candidate spin-3/2 atoms are alkaline-earth atoms of 135 Ba and 137 Ba. Their 
resonances of 6s 2 — > 6s 1 6p 1 are at 553.7 nm^S, thus making them as promising 
candidates for laser cooling and further experimental investigation. Their scattering 
lengths are not available now, but that of 138 Ba (spin 0) was estimated as — 4Ias 
ESI. Because the 6s shell of Ba is fully-filled, both the a , a 2 of 135 Ba and 137 Ba 
should have the similar value. Considering the rapid development in this field, we 
expect that more spin-3/2 systems can be realized experimentally. 

In the rest of this paper, we will review the progress in spin-3/2 system, and also 
present some new results unpublished before. In Sect. [21 we introduce the proof of 
the exact 5*0(5) symmetry. In Sect. [31 we present the effect of the 50(5) symmetry 
to the spin-3/2 Fermi liquid theory, the mean field phase diagram of the lattice 
Hubbard model, and the sign problem of quantum Monte Carlo sign problem. In 
Sect. 21 we review the quintet Cooper pairing and the topological defect of the half- 
quantum vortex. In Sect. Owe discuss the quartetting instability and its competition 
with pairing instability in ID systems. In Sect. [SI we study the spin-3/2 magnetism. 
We summarize the paper in Sect. [7] 

Due to limitation of space, we will n ot c over several important works on the 
spin-3/2 fermions. For example, Hattori S3] studied the Kondo problem finding 
an interesting non-Fermi liquid fixed point. Recently, the ID spin-3/2 model with 
arbitrary interactions was solved using Bethe-ansatz by Controzzi et. a/SH. 



2. The hidden SO(5)(Sp(4)) symmetry 

In thi s section, we will review the hidden symmetry properties of spin-3/2 systems 
21 23 -[jet us first look at a familiar example of hidden symmetry: the hydrogen 
atom. The hydrogen atom has an obvious 50(3) symmetry which can not explain 
the large energy level degeneracy pattern of n 2 . Actually, this degeneracy is not 
accidental. It arises from the 1/r Coulomb potential which gives rise to a hidden 
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conserved quantity, the Runge-Lentz (R-L) vector. This vector describes the orien- 
tation and eccentricity of the classical elliptic orbits. The R-L vector and orbital 
angular momentum together form a closed SO (4) algebra which is responsible for 
the degeneracy pattern of n 2 . Now we can ask: what are the hidden conserved 
quantities in spin-3/2 systems? We will answer this question below. 

2.1. Model Hamiltonians 

We st art with the generic form of the spin-3/2 Hamiltonian of the continuum model 
| 45 | 20 | rp^g om y assumptions we will make are the spin SU (2) symmetry and the s- 
wave scattering interactions. Required by Paul's exclusion principle, the spin chan- 
nel wave function of two fermions has to be antisymmetric in the s-partial wave 
channel. As a result, only two independent interaction channels exist with the total 
spin F = 0, 2. Taking the above facts into account, the Hamiltonian reads 

H= fd d r{ - - ^ a (f) + g Pl (r)P , (r) 

a=±3/2,±l/2 

+ 92 £ Pl m (f)P*,mW}, (1) 
TO=±2,±1,0 

with d the space dimension, pt the chemical potential. P§QiP\ m are the sin- 
glet (St = 0) and quintet (St = 2) pairing operators defined through 
the Clebsh-Gordan coefficient for two indistinguishable particles as Pp m (r) = 

E a Sh F ^ m \U a ^U^U^' where F = °' 2 and m = ~F,-F + 1,...,F. The 
feature that only two interaction channels exist is important for the existence of 
the hidden SO (5) symmetry. 

If spin-3/2 atoms are loaded into optical lattices, we need to construct the lattice 
Hamiltonian. We assume Vb the potential depth , k = tt/Iq the light wavevector, and 
lo the lattice constant. The hopping integral t decreases exponentially with increas- 
ing Vb- Within the harmonic approximation, U/AE ss (n 2 /2)(a a /lo)(Vo/E r ) 1 / 4 , 
with U the repulsion of two fermions on one site, AE the gap between the lowest 
and first excited single particle state in one site, a s the s-wave scattering length in 
the corresponding channel, and E r = h 2 k 2 /2M the recoil energy. In the absence 
of Feshbach resonances, the typical estimation reads a s ~ lOOos (as the Bohr 
radius), l - 5000A, and (V /E r )^ 4 » 1 ~ 2. Thus we arrive at U/AE < 0.1, and 
the system can be approximated by the one-band Hubbard model 

H = ~ l H { C 'W + h - c -} -»I2 c lci<? + UoY, p o,o(i)Po,o(i) 

(ij) ,<r icr i 

+ U 2 P lm(i)P2,m(i)- (2) 

i,m=±2,±l,0 

for the particle density n < 4. At half-filling in a bipartite lattice, fi is given by 
A*o = (Uo + bU2)/4 to ensure the particle-hole (p-h) symmetry. The lattice fermion 
operators and their continuum counterparts are related by ip a (r) — c a (i) / '(lo) d ^ 2 ■ 
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E,=0 

E 4 =Lj,-2n - 
Ef =U ()+ 5^-4n 

Fig. 1. Energy level diagram for a single site problem. The longer and shorter arrows denote the 
spin components of Sz = ±| and S z = ±^ respectively. The sets of Ei,4,6 (spin singlet) are also 
SO (5) singlet; E2,s(spin quartet) are SO(5) spinors; E3 (spin quintet) is an SO(5) vector. 



The best way to understand the SO (5) symmetry without going through the 
mathematical details is to look at the energy level diagram for the single site prob- 
lem as depicted in Fig. [I] It contains 2 4 =16 states which can be classified into three 
sets of spin singlets (-©1,4,6), two sets of spin quartets (©2,5), and one set of spin 
quintet {E 3 ). Interestingly, this energy level degeneracy pattern matches SO (5) rep- 
resentations perfectly: the spin singlets ©1,4,6, quartets -©2,5, and quintet E3 can be 
considered as 50(5) singlets, spinors, and vector as well. Thus the single site spec- 
trum is SO(5) symmetric without any fine tuning. If we further tune Uq = U2 = U, 
then £3 4 become 6-fold degenerate making the system SU(4) invariant. In this case, 
interactions in Eq.[2]reduce to the density-density interaction as Hint = ^n(n— 1). 
The ST/ (4) symmetry simply means that the four spin-components are equivalent 
to each other. 



2.2. The particle-hole channel SO (5) algebra 

In order to prove the 15*0(5) symmetry in the many-body Hamiltonians of Eq.[T]and 
[U we need to construct the SO (5) algebra. The four spin components ip a (f)(a = 
±|, ±i) render 4 2 = 16 p-h channel bi-linears. As a result, the charge operator and 
three spin F X . V:Z operators are not a complete set for describing the total degrees of 
freedom. High order tensors are needed in terms of the tensor product of s = 3/2 
spin matrices Fj: 



rank-0: I, 

rank-1: F\ £=1,2,3, 

rank-2: QFiFj, a £?• = $ = 0, 

rank-3: ^Jifj-Ffc, L=l,..,7, $ fc = $ fc , tf ik = 0, (3) 



I 



l 
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where £?• and £°- fc are fully symmetric, traceless tensors with the detailed forms 
given in Ref. SSI 

The five rank-2 tensors are often denoted as spin-nematic matrices. They are 
T 1 = -j=(F x F v + F y F x ), T 2 = -j=(F,F s + F X F Z ), T 3 = ^={F s F y + F y F z ), 

r 4 = ^(^ 2 -^), T* = ±{l%-F*). (4) 

Remarkably, for F = |, they anticommute with each other and form a basis of 
the Dirac T matrices as {F a ,rfc} = 26 a b, i.e., they form an 5*0(5) vector. More 
explicitly, they are expressed as 

where / and a are the 2x2 unit and Pauli matrices. On the other hand, three spin 
operators F XtVtZ and seven rank-3 spin tensors ^ k FiFjFk,(L — 1,..,7) together 
form the ten 5*0(5) generators. By linear combinations, these generators can be 
organized into 

r ah = -^[r a ,r b ] (i<a,&<5). (6) 

Consequently, the 16 bilinear operators can also be classified according to their 
properties under the 50(5) transformations as scalar n (density), vector n a (spin 
nematics), and anti-symmetric tensors L a b (spin and spin cubic tensors): 

n{r) = ^(r)^(r), n a (f) = -^(^^(r), 

L ab (f) = (7) 

L a b and n a together form the SU(4), or isomorphically, the 50(6) generators. The 
n, n a , and L a b in the lattice model can also be defined accordingly. 

Next we study the particle-particle channel bilinears. Pairing operators can be 
organized as 50(5) scalar and vector operators through the charge conjugation 
matrix R = T1T3 as 

rp{f) = Re?7 + i Imrj = ^^(f)i? Q ^^(r), 

X\ (r) - Re Xa + 1 Im Xa = ^i{r){Y a R) a ^\(r). (8) 

The quintet pairing operators \\ are J us t the polar combinations of F z 's eigenop- 
erators P\ m with the relation 

pt __V ] P t _ TXi + ixj p t _ -X3 ± ixl p t _ ,xl /q\ 
~~ ' ±2 ~~ 2 ' — 2 ' ,0 ~~ v2' 
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The existence of the R matrix is related to the pseudoreality of 50(5) 's spinor 
representation. It satisfies 

R 2 = -1, i? f = iT 1 = l R = -R, RT a R = -*r°, RT ab R = t T ab , (10) 

where t Y a ' ab are the transposed matrices of T a ' ah . The anti-unitary time-reversal 
transformation can be expressed as T = R C, where C denotes complex conjugation 
and T 2 = — 1. N, n a , and L a }, transform differently under the T transformation 

TnT^ 1 = n, Tn a T~ x = n a , TL ah T~ x = -L ab . (11) 

Under particle-hole transformation, fermions transform as %p a — > R a ptl)^,^ a — > 
Rapipp- Correspondingly, L a b,n a transform as 

With the above preparation, the hidden 50(5) symmetry becomes manifest. We 
can explicitly check all the ten SO (5) generators L a (, operators commute with the 
Hamiltonian Eq. [1] and Eq. O In other words, the seven rank-3 spin tensor opera- 
tors are hidden conserved quantities, which play the same role to the Runge-Lentz 
vectors in the Hydrogen atom. The kinetic energy parts in Eq. [1] and Eq. [2] have 
an explicit SU (4) symmetry. The singlet and quintet interactions are proportional 
to rf(r)r](r) and Xa(^)Xa(r) respectively, thus reducing the symmetry group from 
SU(4) to SO(5). When g = .92 or U = U 2 , the SU(4) symmetry is restored be- 
cause XaiV together form its 6 dimensional antisymmetric tensor representation. 
In the continuum model, interactions in other even partial wave channels also keep 
the S*0(5) symmetry. The odd partial wave scattering include spin 1 and 3 channel 
interactions g\ and 53, which together could form the 10-d adjoint representation 
of 50(5) at gi = gz- However, to the leading order, p-wave scattering is weak for 
neutral atoms, and can thus be safely neglected. In the lattice model, this corre- 
sponds to the off-site interactions, and can be also be neglected for neutral atoms. 
The proof of SO(5) invariance in the continuum model applies equally well in the 
lattice model at any lattice topology and at any filling level. For later convenience, 
we rewrite the lattice Hamiltonian in the following form: 

Hj = - J2 {-^(«W-2) 2 + ^( l )}-( M - Mo )^n( ? ), (13) 

2,l<a<5 i 

with V = (3J7 + 5f7 2 )/8 and W = (U 2 - U )/2. 

It is natural to think about the symmetry properties in other high fermionic 
systems with F = n — i (n > 3). Indeed, we can generalize this Sp(4) symmetry to 
the Sp(2n) symmetry under certain conditions, but fine tuning is generally speaking 
needed. In such systems, the s-wave channel interactions can be classified into 
channels with total spin St = 0, 2, .., In — 2 with coupling constants (?0i gi-, <?2ri-2, 
respectively. The Sp(2n) symmetry appears when g 2 — 94 = ■■■ = Qin-i- The proof 
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of the Sp(2n) symmetry and an introduction to the Sp(2n) algebra are given in 
Ref.E3] 



2.3. The SO(8) structure in the bipartite lattice systems 

We next further explore the full symmetry of the spin-3/2 system by extending the 
above algebra to include the p-p channel. In spin-1/2 bipartite lattice systems, it 
is well known tha t the particle density and ry-pairing operators form a pseudospin 
SU(2) algebraS3. For spin-3/2 systems, we embed the above 50(5) algebra inside 
the 50(8) algebra involving both p-h and p-p channels. 50(8) is the largest Lie 
algebra which can be constructed with four component fermions. The 50(8) sym- 
metry is shown to be dynamically generated in spin-i two-leg ladder systems, and 
its effect was extensively studied &8|49 [ <j,-^ e 50(8) algebra here is attributed with 
different physical meaning. Its generators M ab (0 < a < b < 7) reads 



/ Rexi 



M, 



ab 



Rexs N 
Imxi 





Re?y \ 

Til 



V 



n 5 

-Im/y 




with N = (n— 2)/2. In the bipartite lattice, we define the global 50(8) generators as 
M a b = J2i M ab{i) for M ab = n, n a , L ab in the p-h channel, and M ab = Yji(-Y M ab{i) 
for M ab — Re?7, Imry, Rexa, Imx a in the p-p channel. 

At Uq — 5 U 2 and a = u (fi = Uo+ ^ U2 ), the three spin singlet sets of Ei ji>e be- 
come degenerate. They form a spin-1 representation for the SU(2) algebra spanned 
by ,r),N. In this case, Hj can be rewritten as 



Hi= J2 {-U 2 Ll b (i)-(fM-f, )n(i)}, 



(14) 



i, l<a,6<5 



by using the Fierz identity Ei< a <s L ab(*)+Ei<a<5 n 2 a (i)+5N 2 (i) = 5. The symme- 
try at half-filling is 50(5) £§> SU(2), which unifies the charge density wave (CDW) 
and the singlet pairing (SP) order parameters. Away from half-filling, this symme- 
try is broken but rj^rf are still eigen-operators since [H, rf] = — (p — po)rj and 
[H,V\ = if* - Vo)v- 

Similarly, at Uq = —3U2 and p — fiQ, the two spin singlet states E\$ and the 
spin quintet states become degenerate. They form a 7- vector representation for 
the 50(7) group spanned by M a b(0 < a < b < 6). The H] can be reorganized into 

]T {~U 2 M ab (i) 2 -(u-u )n(i)}. (15) 

i,0<a<b<6 



H, = 



The 50(7) symmetry is exact at half-filling. Its 7-d vector representation uni- 
fies the order parameters of the staggered J2i(~ Y n a{i), and the singlet pairing 
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Si *7 (*) • ^ s 21-d adjoint representation unifies the order parameters of the stag- 
gered J2i(~T L ab(i), the CDW SiC - )*-^), and quintet pairing J2iXHi)- Away 
from half-filling, quintet pairing operators are spin 2 quasi-Goldstone operators 
[H, xt(Xa)] — T(m — MqIXo(Xo)- These x- modes are just the analogs of the n modes 
in the high T c context. 

3. Effects of the SO(5) symmetry 

In this section, we review the SO(5) effect to the Fermi liquid theory the mean 
field p hase diagram 1211, the sign problem in quantum Monte-Carlo simulations 
21 | 23 [ -y^g a ^ so construct effective strong coupling models at half- filling in Sect. 
1331 

3.1. The SO(5) Fermi liquid theory 

The Fermi liquid theory in spin 3/2 systems is simplified by the 50(5) symmetry. 
A particle-hole pair in such systems can carry total spins with St = 0, 1, 2, 3, thus 
it generally requires four Fermi liquid functions in these channels. However, the 
5*0(5) symmetry reduces them to three independent sets as 

pa pa pab pao 

f a p,7s(i>,&) = f*(p,fi) + fv{p^){^-) a p{^-) 1 s + ft{p,0)(— )cp{— )jS, (16) 

where p and pi are the momenta of two particles near the Fermi surface, f s<v t 
describe the interaction in the 50(5) scalar, vector, and tensor channels, re- 
spectively. Within the s-wave scattering approximation, these functions become 
constants as f s = (g + 5g 2 )/l6, fv = {go - 3g 2 )/4, ft = -(go + ff2)/4. In 
other words, the two channels of St — 1,3 are degenerate and included in the 
tensor channel. The susceptibility in each channel reads Xs.v,t{r — f*,t — t) = 
-i8(t-t')(\[d SyV ,t(r,t),d s , Vyt (r / ,t')]\) with StV>t = n, n a , L ab , respectively. Within 
the random-phase approximation, the Fourier transforms in the momentum and 
frequency space are Xs,v,t(q,v) = X°(<7,^)/[1 + fs,v,tX°{<l, v)], where X°(g,o;) is 
the standard Lindhard function of the free systems. The degeneracy in the spin 1 
and 3 channels was first pointed out in Ref. 1^*1, but it appears accidentally there. 
It is actually exact and protected by the generic 50(5) symmetry. Experiments 
in the Fermi liquid regime can determine the four Fermi liquid constants in the 
St = 0, 1, 2, 3 channels separately and verify the degeneracy between spin 1 and 3 
channels. 

3.2. Mean field analysis at half-filling in bipartite lattices 

In the weak coupling limit, we perform a mean-field analysis to the lattice Hamil- 
tonian Eq. [5] at half-filling in a bipartite lattice. The decoupling is performed in all 
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Fig. 2. The MF phase diagram at half- filling in a bipartite lattice. A) and B): staggered phases of 
the 50(5) adjoint and vector Reps; C): the singlet superconductivity; D): CDW; E), F), G) and H): 
exact phase boundaries with higher symmetries. Between the dashed lines (Uq < U2 < — 3/5Ub) , 
a Monte-Carlo algorithm free of the sign problem is possible. From Wu et al. Ref. [21]. 



of the direct, exchange, and pairing channels as 

H MF = -t { c L c j> + h - c ] - (M - Mo) n W ~ W2 2 ~ 51 ( n a(«)^a(«) 

(ij,cr) i,& i,l<a<5 

i,l<a<fc<5 i 

+ U J2 {< Re ^)> Rer;(z) + (Imr;(z)) Im^i)} + U 2 £ { (Re XQ (i))Re X a(«) 

£ ?.l<a<5 

+ (ImXa(i))Im X a(i)} (17) 

We solve the Ea . [T"7l self-consistentlv at half-filling in the 2D square lattice with 
the phase diagram depicted in Fig. [2] We use the following MF ansatz: 

(n.(i)) = (-Tn a , (N(i)) = (-YN, (L ab (i)) = {-) l L ab 
(»?(*)> =W, (Xa(*')>=5Ca- (18) 

There are four phases A, B, C and D with bulk area separated by high symmetry 
lines E, F, G and H. In phase A, the staggered order parameters ((— ) % L a b(i) 7^ 0) 
form 50(5)'s adjoint representation (Rep) with the residue symmetry 5*0(3) ® 
50(2), thus the Goldstone (GS) manifold is 50(5)/[50(3) ® 50(2)]. In phase B, 
the staggered order parameters ((— ) l n a (i)) ^ form 50(5)'s vector Rep with 
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Phase 


order parameters 


GS manifold 


GS modes 


A 


{-YL ab {i) 


50(5)/[50(3) ® 50(2)] 


6 


B 


(-)*«-(<) 


50(5)/50(4) = 5 4 


4 


C 


n(i) 


[7(1) 


1 


D 


CDW 


1 


/ 


E 


(-)*n (i), (-)*L o6 (i) 


U(A)/[U{2)®U(2)\ 


8 


F 


(-) l n a (i),v(i) 


50(7)/50(6) = 5 6 


6 


G 


CDW, n(i) 


50(3)/50(2) = S 2 


2 


H 


CDW, xa(i), R%&« 


50(7)/[50(5) x 50(2)] 


10 



Tabic 1. Order parameters, the corresponding Goldstone manifolds and the number of Goldstonc 
modes in each phase on the bipartite lattice at half- filling. From Wu et al. Ref. [21]. 







4 


4 




4- 








M- 


4 


4f 



Fig. 3. Three configurations of classical Neel states on the SU(A) line. 

the residue symmetry 50(4) and the GS manifold 5* 4 . In phase C and D, order 
parameters are the singlet superconductivity (rj(i)} ^ and charge density wave 
(CDW) ((—) l N(i)) 7^ respectively, where the SO (5) symmetry is unbroken. SU(A) 
symmetry arises at line E where the staggered order parameters (— ) 8 L a &(i) and 
(-Yn a (i) become degenerate with GS manifold U(4)/[U(2) <g> [7(2)]. Line F and H 
are characterized by the 50(7) symmetry. On the line F, order parameters (— ) l n a {i) 
and rj(i) are degenerate with the GS manifold 5 6 . One the line E, order parameters 
(-YLabii), CDW and X a {i) are degenerate with the GS manifold 50(7)/[50(5) <g> 
50(2). The 50(5) <g> SU{2) symmetry arises at line G, where CDW and singlet 
pairing become degenerate. Order parameters in each phase and corresponding GS 
modes are summarized in Table 3.1. 



3.3. Strong coupling modes at half-filling 

Next we construct the effect Hamiltonians for the bosonic sector in the strong 
coupling limit. Around the 5[7(4) line E with Uq = U2 = U, the low energy states 
are the six states with double occupancy £^4. The exchange model in this projected 



I 



I 
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Hilbert space reads 

H e ff = JL L ab{i)L ab (j) + J n ^ n a(i)n a (j) + E L 2 ab (i), (19) 

l<a<b<5 l<a<5 i,l<a<b<5 

with A[/i = U 2 - U Q . Along the SU(4) line, AUi = and J L = J n = 4^, Eq. 

fT9l is just the 577(4) Heisenberg model with each site in the 6-d representation. At 

three dimensions, we expect the ground state to be long range Neel ordered. Three 

possible classical configurations A, B, C are shown in Fig. [3l which can be rotated 

into each other under SU (4) operations. They are eigenstates of three SU (4) Cartan 

generators: L9.1 = ■kins — m + n-i — n-3), L15 = 4 (713 + m — n-i — n-3), and 
" * v 2 2 — ~ z ^ 2 2 — ~ 

TI4 = \{ns -Tii - n-i + n-3 ). A nonzero onsite AXJ\ term reduces the symmetry 
down to 50(5) and brings the difference between Jl and J n . The latter contribution 
is at the order of (t/U) 2 compared to the former, thus is less important. In phase A, 
five quintet states become the lowest energy states, then the J n and AUi terms can 
be neglected, and thus Eq. [19] reduces to the SO(5) Heisenberg model. The Neel 
states in phase A can be represented in Fig. [3] A and B which breaks TR symmetry. 
In phase B, the singlet state becomes the lowest energy state, Eq. [3] reduces to 
the 5*0(5) rotor model whose (spin-nematic) Neel order keeps the TR symmetry as 
represented in Fig. [3]C. A transition from a site singlet phase to the Neel ordered 
phase happens at zJ\ ~ AUi where z is the coordination number. 

Next we look at the SO(7) line F where the lowest energy state is the singlet 
state £4, and the next bosonic states are the 7- fold degenerate sets of -Eb,3,6- In 
phase B or C , £0,6 becomes higher or lower than £3, respectively. An anisotropic 
SO (7) rotor model can be constructed to describe the above effect: 

H efft2 = J n Y n a{i)n a {j) - y-jy^Mi) + vi^rfti)} 

l<a<5 



AC/2 E M^) + ^N^), (20) 

0<a<6<6 



2 



where AU 2 = lh ^ L and AU Z = t/o+ 2 3C/2 . Along line F, AU 3 = and J n = J v , thus 
the antiferromagnetic spin-nematic order and singlet superconductivity become de- 
generate. This 5*0(7) symmetry is the generalization of the 50(5) theory of high T c 
superconductivity I n phase C, AC/3 < 0, the staggered spin-nematic order n a is 
disfavored, and Eq . [201 reduces to the U(l) rotor model for the singlet superfluidity. 
On the other hand, for AC/3 > in phase B, Eq. [50] reduces to the 50(5) rotor 
model for n a . 

Along line G and in phases C and D, the lowest energy states are three 50(5) 
singlets of £1,4,6 which form a spin-1 Rep for the pseudospin 5C/(2) algebra. The 
effective model is anisotropic spin-1 Heisenberg-like as 

Heff,3 = ^WiMj) + v(i)vHj)} + JnN(i)N(j) + AC/ 4 /V 2 «, (21) 
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with AUi = - U °+ 5U 2 . Along the 50(5) S77( 2 ) line G, J,, = J„ = ^ and 
AC/4 = 0, the singlet superfluidity and CDW are degenerate. In the phase O where 
AC/4 > 0, Eq . [2T1 reduces into the U(l) rotor model for superfluidity, while in the 
phase D where AU4 < 0, Eq. [21] reduces into the Ising model for CDW. 

Similarly, around the other 5*0(7) line H, the sets of -Ea,4,6 consist the low 
energy Hilbert space. We can construct the anisotropic SO(7) Heisenberg model in 
the 7-vector representation as 

Heff = J L £ L ab {i)L ab (j) - y J2 {xUi)XaU)+Xa(i)xi(J)} 

l<a<b<5 l<a<5 

+ J N N(i)N(j)+AU 5 ( 22 ) 

i,l<a<b<5 

where Jl — J x = Jn = ut+u^ an< ^ AC/5 = around line H. In phase A where 
AC/5 > 0, Eq. [22] reduces into the 50(5) Heisenberg model, while in phase B where 
AC/5 < 0, it reduces into the Ising model for the CDW order. 

3.4. Quantum Monte-Carlo sign problem in spin-3/2 systems 

The sign problem is a major difficulty for the quantum Monte-Carlo (QMC) simu- 
lations to apply to fermionic systems. In spin-3/2 systems, due to t he structure of 
50(5) algebra, the sign problem in the the auxiliary field QMC^Hlis shown to be 
absent in a large part of the phase diagram I 21 | 51 | 23 ] 

We use the equivalent version of the lattice Hamiltonian Eq. [13] where inter- 
actions are represented in terms of the TR even operators n and n a . By using 
the Hubbard-Stratonovich (H-S) transformation, the four-fermion interactions of 
n 2 and n\ are decoupled at V, W > (or —3/5 Uq > U% > Uq), and the resulting 
partition function can be written as 

Z = cxp{-^ f dT^2{n(i,T) - 2) 2 - W f 'dT^n^(i,r)|det|/ + fij. (23) 

The functional determinant I + B = I + Te~ fo dTH K+Hi(T) - g f rom ^h e integration 
of fermion fields, n and n a are real H-S bose fields. The kinetic energy Hk and 
time-dependent decoupled interaction Hi(r) read 

h k = -t^2(i>l a ipj,a + h.c.) - ^(u - m)*l>Lipi<T 

i i 

flint (t) = -Vj2i>iAT)i>iAr)n(i,T) - W^24 !a (T)T%^ il(T >(T)na(i,T). (24) 

i i,a 

Generally speaking, det(/ + B) may not be positive definite, making it difficult 
to use the probability interpretation of the functional integrand in the QMC sim- 
ulation. However, if the above decoupling scheme satisfies an TR-like symmetry, 
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the following theorem states the positivity of det(7 + B) I 52 | 53 | 23 [ ^ th ere exists an 
anti-unitary operator T, such that 

THkT- 1 = H K , THiT~ x = Hi, T 2 = -1, (25) 

then the eigenvalues of the I + B matrix always appear in complex conjugate pairs 
no matter I + B is diagonalizablc or not, i.e., if A^ is an eigenvalue, then A* is 
also an eigenvalue. If A^ is real, it is two-fold degenerate. In this case, the fermion 
determinant is positive definite, 

dct(/ + £) =]J|A. 1 | 2 >0. (26) 

i 

The detailed proof can be found in Ref. The T operator does not need to be the 
physical TR operator. As long as Eq. [2Uis satisfied, the positivity of det(7 + B) is 
guaranteed. We emphasis that the above criterion applies for any lattice geometry 
and doping levels. 

Because both the density n and spin-nematic operators n a are TR even oper- 
ators, the sign problem disappears as long as V, W > or —3/5 Uq > Ui > Uo 
(see Fig. [5J . This region includes interesting competing orders such as the staggered 
spin-nematic order n a , singlet superfluidity 77, CDW, and also the phase boundaries 
of the 50(7) line F and the SO{5) ® £17(2) line G. The absence of the sign prob- 
lem provide a good opportunity to study the competition among these orders, in 
particular, then doping effect, the frustration on the triangular lattice, etc, which 
are difficult at low temperatures for previous Monte-Carlo works. 

The above sign problem criterion can also be applied to other systems, such as 
multi-band Hubbard models. For example, Capponi et al. showed a ground state 
staggered current order b y QM C simulations without the sign problem in a bi-layer 
extended Hubbard model ^4 



4. Quintet Cooper pairing and half-quantum vortices 

When the interaction constant in the quintet channel 92 goes negative, spin 3/2 
systems can support quintet Cooper pairing, i.e., Cooper pairs with total spin 2. 
In this section, we discuss the collective modes in this state, and the topological 
defect of half-quantum vortex and related non-Abelian Cheshire charges L ^- ±l . 

4.1. Collective modes in the quintet pairing states 

We write the BCS mean field Hamiltonian for the quintet Cooper pairing state as 

^2v72 



Hmf = f d D r{ ]T (r) (_|J_ _ a ) M r) + £ xUr)A a (r) + h.c. 
- iA;(r)A,(r)}. (27) 

09 J 
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A a is proportional to the ground state expectation value of the quintet pairing 
operators \a by A a (r) = <72(xl( r )) (a = 1 ~ 5). The five \a operators de- 
fined in Eq. [8] are the spin channel counterparts of the five atomic d-orbitals 
(d xy , d xz , d yz , d^ z 2_ r 2, d x i_ y i). They transform as a 5- vector under the 50(5) group. 

The Ginzb urg- Landau (GL) free energy around T c for the quintet pairing was 
studied in Ref. without noticing the hidden 50(5) symmetry. Facilitated by this 
high symmetry, we organize the GL free energy in an explicitly 50(5) invariant 
way, and give a general analysis to the quintet pairing structures. The GL free 
energy reads 

F GL = Jd D r 7 VA:VA Q + a(r)A;A a + ||A:A a | 2 + | £ | KU ^ 1 ~ W 

ldn. T Q a ldn7C(3) nh 2 7((3) , OD , 

a = ~2Te {l -T} ^ 1 = ^ 2= 2&8^f' 7= 4M8^f (28) 

where n is the particle density, dn/de is the density of states at the Fermi level, 
and ((x) is the Riemann zeta function. At /3 2 > or /3 2 < 0, two different pairing 
states are favorable to minimize the free energy, i.e., the unitary polar state and 
non-unitary ferromagnetic state respectively. The polar state is parameterized as 
A a = \A(T)\e l6 d a , where 9 is the U{1) phase, and d = d a e a is a 5-d unit vector 
in the internal spin space. The ferromagnetic pairing state is described as A a = 
|A(T)|e* e (di a + i d 2a ) with d\, d 2 two orthogonal 5-d unit vectors. The standard 
Gor'kov expansion gives 02 > 0, thus the polar state is stable in the framework of 
the BCS theory. In the following, we focus on the polar state which is also expected 
to be stable at zero temperature. 

At T = OK, various low energy dissipationless Goldstone modes exist due to 
the broken symmetries. In addition to the usual phonon mode, four branches of 
spin wave modes carrying S = 2 arise from the breaking of the SO(5) symmetry 
down to 50(4). The above GL free energy constructed around T c does not apply 
to describe these Goldstone modes. For this purpose, Cooper pairs can be treated 
as composite bosons. This treatment gives a g ood de scription of the phonon mode 
in the neutral singlet BCS superfluid in Ref. I 54 | 55 | jj ere) we generalize it to the 
quintet pairing by considering a phenomcnological Hamiltonian for spin 2 bosons 

Heff =Jd D r^Y, I W °! 2 + 2£-(*£*» - + ^~ E(**^*M 2 , (29) 

a " S P a<b 

where ^ a s are the boson operators in the polar basis, — i(S ac 5bd — $addbc) are the 
50(5) generators in the 5x5 representation, the equilibrium Cooper pair density po 
is half of the particle density Pf,Xp an d Xsp are proportional to the compressibility 
and 50(5) spin susceptibility respectively. Taking into account the Fermi liquid 
correction, x P = iV//[4(l + Fo)] and Xs P = 2V//[4(1 +i*o)], where Nf is the fermion 
density of states at the Fermi energ y, lJT o' are the Landau parameters defined in 
the 50(5) scalar and tensor channels ^ respectively. We introduce p(r) and l a b( r ) 
as the Cooper pair density and 50(5) spin density respectively, and parameterize 
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Fig. 4. A) The Goldstone manifold of d is a 5D hemisphere RP 4 . It contains a class of non- 
contractible loops as marked by the solid curve. B) The 7r-disclination of d as a HQV. Assume 
that d || at </> = 0. As the azimuthal angle <j> goes from to 2ir, d is rotated at the angle of rf>/2 
around any axis h in the 5 3 equator spanned by ei^ 3,5- From Wu et al. Ref. [24], 



ty a = y/~p e d a - Using the standard commutation rules between p and 9, l ao and 
d a , we arrive at 

x P d?e - = o, (so) 

h 2 p s 

d t lab = -^-^-{da^db- d b W 2 d a ), Xspd t d a = -l ab d b , (31) 

where p s and p sp are the superfluid density and spin superfluid density respectively. 
At T = OK in a Galilean invariant system, p s is just P t/% while p sp receives Fermi 
liquid corrections as p sp /p s = (1 + F" /3)/ (l + F*/3) ^ where is the Landau 
parameter in the SO(5) vector channel ^J. The sound and spin wave velocities are 
obtained as v s = -J po/(2x p M) and v sp = y i /po/(2x sp M), respectively. 



4.2. Half-quantum vortex 

Superfluids with internal structure can support half-quantum vortex (HQV) as a 
stable topological defect because the order parameter contains a Z 2 gauge symme- 
try. For the quintet pairing case, A a = |A|e l6l a! a is invariant under the combined 
operation 

d -^-d, 9^8 + 7T. (32) 

As a result, d is actually a directionless director instead of a true vector. Thus the 
fundamental group of the GS manifold is iti{RP a ® U{\)) = Z ® Z 2 as depicted in 
Fig. 0] A. The Z 2 feature gives rise to the existence of HQV as depicted in Fig. 2] 
B. As we move along a loop enclosing HQV, the 7r phase mismatch in the 9 field is 
offset by a 7r-disclination in the d- field, thus A 's are still single- valued. 
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Next we check the energetics of HQV. The static energy function can be written 

as 

E = J d D r^{ Ps (V6f + p sp (Vd) 2 }. (33) 

The energy density per unit length of a single quantum vortex is E\ = p s log — , 
while that of two isolated HQVs is E 2 = ^g(Ps + p sp ) log Thus for p sp < p s , a 
single quantum vortex is energetically less favorable than a pair of HQVs. Although 
at the bare level p s — p sp , p sp receives considerable Fermi liquid correction and 
strong reduction due to quantum fluctuations in the 5D internal space. Generally 
speaking, p sp < p s holds in terms of their renormalized values. ^^^^ 

An interesting property of the HQV is its Alice string behavior I57 | 58 [ j n this 
context, it means that a quasi-particle change its spin quantum numbers after it 
adiabatically encircles the HQV. For example, in the 3 He-A phase, a quasi-particle 
with spin | changes to J, up to a U(l) Berry phase. The HQV in the quintet pairing 
case behaves as a non-Abelian generalization of the above effect. Without loss of 
generality, we assume that d || £4 at the azimuthal angle <f> = 0. As <f> changes from 
to 2ir, d is rotated at the angle of <p/2 in the plane spanned by 64 and n, where h 
is a unit vector perpendicular to 64, i.e., a vector located in the S 3 sphere spanned 
by ei.2,3.5. We define such a rotation operation as U(n,4>/2). When U acts on an 
5*0(5) spinor, it takes the form of U(n,<p/2) — exp{— i\ nh \ }, when U acts on 
an 50(5) vector, it behaves as U{n,cj)/2) = exp{— i&ribL bA } where L ab, s are the 
SO(5) generators in the 5x5 vector representation. The resulting configuration of 
d is d(h,(f)) = U(n,(f>/2)d(h,0) — (cos ie^ — sin^n). As fermionic quasi-particles 
circumscribe around the vortex line adiabatically, at <p — 2n fermions with S z = i § 
are rotated into S z = ±| and vice versa. For convenience, we change the basis 'J/ 
for the fermion wavefunction to (||), | — |), | — \)) T ■ After taking into account 
the 7r phase of the superfluid vortex, ^ transforms by 

* Q -> % - iU(h, ^) a ^ p =(^ W \ * (34) 
where W is an SU{2) phase depending on the direction of h on the S 3 sphere as 

W(h) = ( n3 + in2 - ni -. in5 ). (35) 
\nx — ins n 3 -in 2 J 

The non-conservation of spin in this adiabatic process is not surprising because the 
SO(5) symmetry is completely broken in the configuration depicted in Fig.[3]B. 

4.3. SO(4) Cheshire charge 

Another interesting concept related to the Alice string in the gauge theory is the 
Cheshire charge, which means a HQV pair or loop can carry spin quantum number. 
When a quasi-particle carrying spin penetrate the HQV loop or pair, the spin 
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Fig. 5. The configuration of a 7r-disclination pair 
azimuthal angles and d{r) \\ {=4 as r — > 00. After a 
with S z = ±§ change to S z = ±^ and vice versa 
Wu et at Ref. [24]. 



or loop described by Eq. 1371 (pi 2 an d A</> are 
fermion passes the HQV loop, the components 
with an SU (2) matrix defined in Eq. 1341 From 



conservation is maintain by exciting the Cheshire charge of the HQV loop or pair. 
Below we will construct the SO (4) Cheshire charge in the quintet Cooper pairing 
state. 

We begin with a uniform ground state where d is parallel to 64 axis where the 
50(4) symmetry generated by T a b(a, b — 1,2, 3, 5) is preserved. This 50(4) algebra 
can be reorganized into two inter-commutable SU(2) algebras as 

r x (T[) = i(±r 35 - r 12 ), t 2 (t^) = ~(±r 31 - r 25 ), 

T 3 (T 3 ') = i(±r 2 3-r 15 ). (36) 

7i ; 2.3 an d 2j 2 3 ac ^ i n the subspaces spanned by | ± |) and | ± i), respectively. 
50(4) representations are denoted by \T, T3; T', T3), i.e., the direct-product of rep- 
resentations of two SU(2) groups. 

Now we construct a HQV loop or pair as depicted in Fig. [5] where <f>\_2 are 
azimuthal angles respect to the vortex and anti-vortex cores respectively. The solu- 
tion of the configuration of the d vector is described by the difference between two 
azimuthal angles A0 — 4>2 — <Pi as 

d[n, A<p) — cos _ ^ _e 4 — sm ~2~ ' ' ' 

where h again is a unit vector on the 5 3 equator. This configuration is called a 
phase-sharp state denoted as \n) vt . Because the above 50(4) symmetry is only 
broken within a small region around the HQV loop, quantum fluctuations of n 
dynamically restore the 50(4) symmetry as described by the Hamiltonian 

H ro t= Yl "2?> M ab = i{n a d nb -n b d fla ), (38) 
0,6=1,2,3,5 
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with the moment of inertial I = \sp J dPr p sin 2 Thus the zero modes \n) vort 
are quantized into the global 5*0(4) Cheshire charge state, which is a non-Abelian 
generalization of the £7(1) Cheshire charge in the 3 He- A phase ^^3. The Cheshire 
charge density is localized around th e vorte x loop. In contrast, the Cheshire charge 
in the gauge theory is non- localized tSZEHl. The HQV loop in the SO(4) Cheshire 
charge eigenstates are defined as \TT^;T'T^) vt = J AeS 3 dn Ftt 3 ;T't' (n) \h) v t, where 
Ftt^t't^? 1 ) are the S 3 sphere harmonic functions. 

When a quasiparticle penetrates a HQV loop or pair in the quintet Cooper 
pairing state, quantum entanglement is generate between them in the final state. 
We demonstrate this through a concrete example, with the initial state \i) made 
from a zero charged HQV loop and a quasiparticle with S z — | as 



/ = j dn \n) vt (g) (u c\ + v c_3 (39) 

JneS 3 2 2 

where is the vacuum for Bogoliubov particles. For each phase-sharp state 

\h) v t, the particle changes spin according to Eq. [34] in the final state |/). The 
superposition of non-Abelian phase gives 

|/)= / dn {u {Wl x c\ +W^cl L )+v (W£ci +W? 2 c_i)\\n) vt ®\n) qp 



nes 3 



dn(h 3 — in 2 )\n) vt ®(uc[|u c_i)\il) qp 



dh(fii — ih§)\h) v t ® (u <J_± — v ci )\il) qp . (40) 

AGS 3 2 2 



In terms of the 50(4) quantum numbers, \i) is in a product state as |00; 00)ut 

I 2 2' ^0)gp, 



"-OO)^, and |/) is 



,11 1-1, , 11, ,11H V , 1-1, 
I 2 2 ; 2~Y )vt ® 1 5 2 2 )qp " ! 2 2 ! 2 2 )vt ® 1 ; 2 ~ )qp ' ( } 

In the channel of the (T',Tg), the final state is exactly an entangled Einstein- 
Podolsky- Rosen (EPR) pair made up from the HQV loop and the quasi-particle. 



5. Quartetting instability in ID spin-3/2 systems 

Quartetting instability is a four-fermion counterpart of the Cooper pairing 
in spin-3/2 systems, whose order parameter can be written as O qrt (r) = 
tp\ {r)ij}\ (r)ip t_i (r)ip L 3 (r). The quartet is an SU(4) singlet which can be consid- 
ered as a four-body EPR state. This kind of order is very difficult to analyze 
in high dimensions because O qr t contains four fermion operators. It lacks of a 
Bardeen-Cooper-Schrieffer (BCS) type well-controlled mean field theory. Neverthe- 
less, considerable progress has been made in ID systems. By using Bethe ansatz, 
Scholttmann showed that the ground state of the SU(2N) model is character- 
ized by the formation of the 2N-particle baryon-like bound state. In contrast, the 
Cooper pairing can not exist because two particles can not form an SU(2N) singlet. 
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1: SU(4) A: Luttinger liquid 

Fig. 6. RG flows in the parameter space (g v ,gt) in the spin channel. Combined with K c in 
the charge channel, they determine various phases as: A) the gapless Luttinger liquid phase, B) 
quartetting phase with QLRO superfluidity at K c > 2 or 2k f CDW at K c < 2. C) singlet pairing 
phase with QLRO superfluidity at K c < 1/2 or 4k f CDW at K c > 1/2. They are controlled by the 
fixed points of (0,0), (+oo, +oo) (line 4), and (— oo,+oo) (line 5) respectively. Phase boundaries 
(line 1, 2, 3) are marked with dashed lines. From Wu Ref. [28]. 

The strong quantum fluctuations in the spin channel in ID suppress the Cooper 
pairing. 

In this section we will review the quartetting instability in ID spin-3/2 systems 
including the SU (4) symmetric model as a special case v^ e w jji gnow that 

both quartetting and singlet pairing are allowed in different parameter regimes, and 
study their competitions. 

5.1. Renormalization group and Bosonization analysis 

In ID systems, we linearize the spectra around the Fermi wavevector kj, and 
then decompose the fermion operators into left and right moving parts as 
i>a = ^R,aG %ksX + ^L,aG~ %ksX ■ The right moving currents are also classified into 
SO(5Ys scalar, vector and tensor currents as Jb{ z ) — Q ( z )'0i?,Q( 2; )i Jr{ z ) — 

^rJz)VrAz) (1 < o < 5), Jf{z) = KaW^W (1 < a < b < 5), 
and the left moving currents can be defined accordingly. The low energy effective 
Hamiltonian density TL = TLo + TLint + 'H'mt 1S written as 

Ho = Vf{jJ R J R + \{J a R J a R + J a R b J a R b ) + (R^L)}, 

Hint = -^JrJl + 9vJrJl + QtJjlJ'Li 

H' int = ^JrJr + 9 Xj R J a R + 9 \jtJ a R +(R^ L). (42) 

The Umklapp term is absent in this continuum model, whose effects become im- 
portant in the lattice model at commensurate fillings, and will be discussed in next 
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Fig. 7. Phase diagram in terms of the singlet and quintet channel interaction parameters go and 
<?2- Configurations at quarter filling after the charge gap opens are shown. Phase A is the SU(A) 
gaplcss spin liquid. Each of phase B and C splits to two parts. With g u — > +co, B.l) quartets 
with both bond and charge orders, C.l) dimerization of spin Peierls order. With g u — » — oo, CDW 
phases of B.2) quartets and C.2) singlet Cooper pairs. Boundaries among phases A, B and C are 
marked with solid lines, and those between phases B.l and B.2, and C.l and C.2 are sketched 
with dashed lines. From Wu Ref. [28]. 



section. At the tree level, these dimensionless coupling constants are expressed in 
terms of go, g 2 defined above as 

, .90 + 5^2 , .9o-3g 2 , .9o + .92 
9c = g c = ^ ' 9v = 9 v = ^ ' 9t=9 t = ^ — , (43) 

which are renormalized significantly under the RG process. The chiral terms in H' int 
only renormalize the Fermi velocities, which can be neglected at the one loop level. 
At g v = gt and g' v = g' t , or in other words go = g 2 , the SU(4) symmetry is restored. 
The above Hamiltonian Eq. [42] can be bosonized through the identity 

3 i 

ip a R,L{x) = r] a /V2Traexp{±iy/n((f> a (x) ± 9 a (x))} (a = +-,+-) (44) 

where the Klein factors rj a are Majorana fermions to ensure the anti-commutation 
relation among fermions of different species. Boson fields cf) a and their dual 
fields 9 a are conveniently reorganized into (f> c (6 C ) in the charge channel, and 
4> V {9 V ), <j>tx(dti), 0*2(0*2) in the spin channels via <j) c . v = 5(^3 ± <f>i ± 0_i + 0_|), 
4>ti.t2 — 5(0^ + <f>! + 0-i — 0-i)- Similar expressions hold for 6*'s. In terms 
of these boson fields, the quadratic part of the Hamiltonian density is standard 

(v = C,V,tl,t 2 ) 



n ° = y E { K v(d x e„) 2 + ^(d x ^) 2 } with 



2-Kv f +gl - g v g' g v 

Kv = i ^+g> v + g^ ^V (V ' + 2^ ~ ( 2. )2 ' (45) 
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where 5*1,42 = <?t, g[x ti = 9t- We need to bear in mind that the expression for K„, v v 
can not be taken seriously at intermediate and strong couplings where finite but 
significant renormalizations to K v and v v take place. The non-quadratic terms are 
summarized as 

Hint = 2(^)2 { COS ^r^*! + C0S V^T(f>t2 1 1 {gt + 9v) COS V47T0„ 

+ (gt - 9v) cos \/4n0v \ ~ / 9 * „ cos\/47r0 t i cos\/47r0 t 2, (46) 
J 2{ira) z 

with the convention of Klein factors as 77377177-177-3 = 1. Eq. 1461 contains cosine 

2 '2 — ~2~ \ — 

terms of both <p v and its dual field 9 V . Their competition leads to two different spin 
gap phases with either cos \f^K(j) v or cos \/~Ak9 v pinned as shown below. The charge 
channel c remains quadratic. 

The renormalization group (RG) equation of g c .v.r. can be calculated from the 
standard operator product expansion technique as 

dg v _ 4 dg t _ 1 2 2 dg c _ 

d\n(L/a) ~ 2ir 9v9u d\n(L/a) ~ 2vr ^ 9t + 9vh dln(L/o) ~ ' 1 j 

where L is the length scale and a is the short distance cutoff. Due to the absence 
of the Umklapp term, the charge part g c remains unrenormalized at the one-loop 
level. The ST/ (4) symmetry is preserved in the RG process along the line g v = g t . 
The RG equations can be integrated as \gf — g%\ = c\g v \ 3 / 2 (c: constant) with the 
RG flows as shown in Fig According to the relation Eq. 031 the corresponding 
boundaries are also shown in Fig. [7l 

The phase diagram Fig [6] (Fig. [Z]) contains three phases at incommensurate fill- 
ings. Phase A is the gapless Luttinger liquid phase lying in the repulsive interaction 
regime with g2 < go- All the cosine terms are marginally irrelevant. The leading 
order divergent susceptibility are the 2fcf-CDW order 02k f .cdw — fpRa^Lai an d the 
2&/-SDW orders in the SO (5) vector channel N a — ^ Ra (^ 2 -) a f3'4'Lp, and those in 
the SO (5) tensor channel N ab — ip\ ia (^-)aP' l pL/3- All of them are with the scaling 
dimension of (K c + 3)/4. 

Phase B is characterized by the formation of quartets. This phase is controlled 
by the marginally relevant fixed point g v = g s — > +oo (line 4 in Fig. [5]), and lies 
in the regime where attractive interactions dominates (Fig. [7]) . The spin gap opens 
with pinned boson fields of <f> v , <pti and <pt2- The pinned values can be chosen as 
(4>v) = (<M) = {4>t2) = up to a gauge degree of freedom. In phase B, every 
four fermions first form quartets, then quartets undergo either superfluidity or 
CDW instabilities. Because the average distance between two nearest quartets is 
d = n/kf, the CDW of quartet is of the 2fc/-type N. Their order parameters 
reduce to Oik s ,cdw °c e l ^^ c and O qr t oc e 2l ^ Sc respectively. Checking their scaling 
dimensions Ajv = jK c and A qrt = 1/K C , O qr t wins over N at K c > 2. Previous 
Bethe-ansatz results for the SU (4) case show the quartets formation in Fig. [B]^l 
Here, we extend the quartetting regime to the whole phase B. On the other hand, 
all the pairing operators 77^ and x a ^ decay exponentially in phase B. 
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Phase C is characterized by the formation of the singlet Cooper pairs. The 
parameter region for this phase is symmetric to phase B in Fig.[6]under the reflection 
g v — > —g v . The corresponding fixed point is — g v — g t — > +oo. Different from phase 
B, the dual field 9 V instead of <p v is pinned. Similarly, we can pin the values of 
bosonic fields (4>ti) = (4>t%) = (&v) = to minimize the ground state energy. The 
two leading competing orders in this phase are the singlet pairing and the CDW of 
pairs. This CDW is just the Akf CDW defined as O^kf^dw = ^Ra^Rp^LpipLa- The 
expressions of these two orders reduce into oc e" 1 ^ 9 ", and O^k^cdw oc e^^. 
At K c > 1/2, the pairing instability wins over the 4fc/-CDW. 

5.2. Ising transition between quartetting and pairing phases 

Next we study the competition between the quartetting (phase B) and singlet 
Cooper pairing (phase C), whic h c an be mapped into the phase locking problem 
of two-component superfluidity One component is defined as A\ — ipt^s, 

and the other one as A3, = ijA ip -i ■ Then the singlet pairing operator and the 

2 ~~2~ 

quartetting operators are represented as 

7? t = i(A| - A\) oc e 4 ^ 9 " cos V^r, 0\ rt = A\a\ = e 1 ^^ cos 2^4> v , (48) 

where the charge channel boson field 6 C is the overall phase, and the spin channel 
boson field 8 V is the relative phase between two components. While depends on 
the relative phase 6 V , 0^ qrt depends the dual field 4> v . 

The overall phase 9 C is always power-law fluctuating in ID, and does not play 
a role in the transition. It is the relative phase v (<p v ) that controls the transition. 
The effective Hamiltonian for this transition is a sine-Gordon theory containing 
cosine terms of both 9 V and <pv as 

H res — — -^7— „ (Ai cos V&tffiv + A2 cos V^tt9 v ) with 
2{na)' ! 

Ai, 2 = (gt±g v ) ((cos V^(f>ti) + (cos a/4tt '<t>t2))- (49) 

In phase C where Ai > A2, then the relative phase <p v is locked giving rise to the 
pairing order. In contrast, in phase B where Ai < A2, the dual field 9 V is locked 
giving rise to the quartetting order. Along the critical phase boundary line 3 in 
Fig. El Ai = A2 or (g v = 0), none of 9 V and <f> v is pinned. This transition can be 
mapped into a theory of two Ising fields one of which is at the critical point. 
Equivalently, it can be regarded as a model with two Majorana fermions, one of 
which is massless while the other is massive. The Ising order and disorder operators 
are given by cos ^f"K^> v and cos ^pn9 v respectively, both of which have the scaling 
dimension of |. 

This Ising symmetry breaking effect can also be understood as follows. In ad- 
dition to the 50(5) symmetry, spin-3/2 systems have another symmetry in the 
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spin channel U, 



e m 4 7r unc [ er w hich ip a , ry , and O qr t transform as 




(50) 



This Z2 transformation shifts the relative phase 6 V as \/tt9 v — ► y 7 ^" ^ 77 but leaves 
\Ztt(/) v unchanged. Thus this Z 2 symmetry is broken in the singlet pairing phase B, 
but is kept in the quartetting phase C. 

5.3. Discussion of quartetting in high dimensions 

The study of quartetting phase, or the more general multi-particle clustering insta- 
bilities in high dimensions is a challenging problem. Ropke et al. compared the 
instabilities of deuteron (pairing) and a-particle (quartetting) channel in nuclear 
physics using diagrammatic method. They found that at low density or strong 
interaction strength, the a-particle instability wins over the deuteron instability. 
Stepnaneko et al. ^1 constructed a trial wavefunction in 3D to describe the quar- 
tetting order. Lee studied the 2D ground state binding energy of the iV-particle 
cluster in the SU(N) case. 

In high dimensions, the competition between quartetting and pairing superflu- 
idities is more complicated than that in ID. First, the pairing superfluidity breaking 
spin rotational symmetry can not be stabilized at ID, but can exist at D > 2. Sec- 
ond, phase transitions from weak to strong coupling regimes can take place at 
D > 2 unlike in ID interactions always renormalize into strong coupling regime. 
Taking into account these facts, we discuss the symmetry class of these competi- 
tions. We begin with the SU(4) line with go = g% = g < 0. At weak coupling, 
the sextet pairing state (rf and x are degenerate) is stabilized breaking the £7(1) 
charge symmetry and also the SU(A), or isomorphically SO(6) symmetry in the 
spin channel. The residue symmetry is [50(5) k Z 2 ] ® Z 2 . The first Z 2 is a com- 
bined spin-phase operation which flips the direction of the pairing operator in the 
spin channel and simultaneously shift the phase by tt, which is of the same nature 
responsible for the HQV in the quintet superfluid in Sect.|H The second Z 2 is purely 
the phase operation ip a — > —ip a which leave pairing operators unchanged. On the 
other hand, the quartetting state in the strong coupling limit is SU(4) invariant, 
and also breaks the U(l) symmetry. Its residue symmetry is SO(6) ® Z4 where Z4 
operation means ip a — » e lm7T / 2 ip a (m = 1 ~ 4). Thus this phase transition is describe 
by the coset SO (6) ® Z 4 /[SO(5) x Z 2 ® Z 2 ] = S 5 . In other words, the pairing to 
quartetting transition can be viewed as an order to disorder transition of pairing 
operators for their spin configuration on the S 5 sphere. If the spin configuration 
becomes disordered, due to the combined spin-phase Z 2 structure, the phase of the 
pair operators is only well defined modulo of it, thus it means the quartetting order 
appears. If g 2 < go < or go < 92 < 0, then in the weak coupling limit, the quintet 
pairing state , or the singlet pairing dominates. A similar reasoning shows 
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that the coset for the transition from the quartetting state to the quintet or singlet 
Cooper pairing states is S 4 or Z%, respectively. 

6. Magnetic properties in spin-3/2 systems 

In the strong repulsive interaction regime Uq, Ui t and at commensurate fillings, 
the low energy physics of spin-3/2 systems are described by the magnetic ex change 
models. In this section, we review the magnetic properties in such systems^ 23, 
and also present new results in Sect. [6~3"l and[ 



6.1. Magnetic exchange at quarter-filling 

We first construct the exchange Hamiltonian at quarter-filling i.e., one particle per 
site. The total spin of two sites in a bond can be Stot = 0, 1, 2, 3. Using the projection 
perturbation theory, we find the exchange exchange energy in each channel as 

At 2 At 2 

Jo = Tr , J2 = TT> -h = -h = 0, (51) 

where the degeneracy of is a consequence of the SO(5) symmetry. This model 
can be represented in the standard Heisenberg type by using bi-linear, bi-quadratic, 
and bi-cubic terms as H ex = J2ij a 0% ' Sj) + b (Si ■ Sj) 2 + c (Si ■ Sj) 3 , with 
a = — jjg(31Jo + 23J2), b = j^(5Jo + I7J2) and c = j$(Jo + Ji)- More elegantly, it 
can be represented in the explicit Sp(A) symmetric form as 

Hex = ^ { ° 4 2 L ab(i)L a b(j) H ^ — -n a (i)n 

a 

(i)}- (52) 

ij 

Eq. [55] satisfies two different SU(A) symmetries at Jo = Ji and J 2 = 0, respec- 
tively. The first one denoted as SU(4)a below is obvious at Jo = J2 where Eq. 
[52] reduces into the SU(A) Heisenberg model of H su ^ t A = J2ij j^L a b(i)L a b(j) + 

n a (i)n a (j)^- Each site is in the fundamental representation. The second SU(A) 
symmetry exists in the bipartite lattice at J2 = denoted as SU(A)b below. We 
perform a particle-hole transformation in the odd sublattice L' ab = L a t, n' a = —n a , 
i.e., the odd sites are transformed into the anti- fundamental representation. Eg. 1521 
is again SU(A) invariant as 

Hsu^b = E ^{L' ab (i)L ab (j) + n' a (i)n a (j)}, (53) 

but with fundamental and anti-fundamental representations alternatively in even 
and odd sites. This SU(4)b model and the more generally staggered S U(N) H eisen- 
berg model were investigated extensively by using the large N method P4 | 65 | These 
two SU (4) symmetries have dramatically different physical properties. At least four 
sites are needed to form an SU(A) singlet for the SU(4)a, while for SU(4)b, two 
sites are enough to form a bond singlet. 
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6.2. ID phases at quarter-filling 

The phase diagram of Eq. [521 at ID can be obtained by using the bosonization 
method to the Hubbard model. At quarter-filling, the 8k f term Umklapp term is 
important and only involves the charge sector. It can be bosonized as 



where g u is at the order of 0(<7o,sf) at the bare level. At quarter- filling, the spin 
channels behave the same as in the incommensurate case discussed in Sect. 15.11 
The charge channel opens up the gap at K c < 1/2 where g u is renormalized to 
infinity. As a result, the degeneracy between the real and imaginary parts of the 
02k f ,cdw and 04k f ,cdw, are lifted. Their real parts describe the usual CDW orders, 
while their imaginary parts mean the bond orders, i.e., the 2k t and 4fc/ spin Peierls 
orders. 

With the opening of charge gap at quarter-filling, various insulating phases 
appear as given in Fig. Phase A and CI corresponds to the regime described 
by the exchange model Eq. EH with J2 > Jq and J2 < Jo, respectively. Phase A 
( J2 > Jo > 0), including the SU (A) a line, is a gapless spin liquid phase. The SU (4) a 
line was solved by Sutherland ^1 In phase C.l (J2 < Jo), the spin gap opens with 
the developing of the Akf spin-Peierls order, i. e., the dimer order. The transition 
between these two phases is Kosterlitz-Thouless like. Insulating phases involving 
CDW order also exist in phases B.l, B.2 and C.2. The quartetting phase B splits 
into two parts B.l and B.2. In phase B.l, the 5/7(4) singlet quartets exhibit both 
the charge and spin Peierls orders. In phase B.2, the 2k f CDW of quartets becomes 
long range ordered. Similarly, the CDW of singlet pairs becomes long range ordered 
in phase C.2. The boundaries between phase B.l and B.2, phase C.l and C.2 are 
determined by the bare value of g u o = as sketched in Fig. [JJ However, due to 
the non- universal relations between g U Q and go.2, the exact boundaries are hard to 
determine. 

6.3. Discussion of the 2D phase diagram at quarter-filling 

The 2D physics of Eq. [52] is a challenging problem. Nevertheless, good understand- 
ing has been achieved for the SU(4)b line (J2 = 0) in the square lattice where 
both quantum M onte-C arlo simulations and the large- TV analysis show that 
the ground state I 65 | 68 | j g \ on g ran g e Neel ordered. As a result of strong quan- 
tum fluctuations, the Neel moments are (— ) % n± = (— ) l Li 5 — (— ) I L 34 w 0.05, 
which is really tiny compared to the SU(2) case. The Goldstone manifold is 
CP(3) = U(A)/[U(1) ® U(3)\ with 6 branches of spin- waves. 

However, the physics is less clear away from the SU (4)b line. The SU(4)a line of 
Jo = J2 is of particular interest. In order to form a singlet for SU (4) a, we need four 
sites around a plaquette with the wavefunction of ^7 e a ^ 7< 5Vi(l)V'^(2)^(3)'(/'|(4)|ri) 
where |f2) is the vacuum. Indeed, an exact diagonalization by Bossche et al. in a 4x4 
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Sl^(4) 
Jo 




Fig. 8. Speculated phase diagram at quarter-filling in a 2D square lattice. Phase A is the plaquette 
ordered phase which is suggested along the SU(A)a line in an exact diagonalization work by 
Bossche et al. 67 Phase C is the Neel ordered phase which is confirmed along the SU(4)b line in 
QMC simulations by Harada et al. 65 We speculate a dimerized phase B lying between phase A 
and C. Dimers in phase B carry antiferromagnetic spin-nematic order. 



sites lattice ^ suggests that the ground state along the SU(4)a line exhibits such 



plaquette order as depicted in Fig. [5] This plaquette order was also found in a large- 
TV mean field theory by Mishra et al<^. However, due to the sign problem in QMC, 
a large size simulation is difficult to confirm this result. Recently, Chen et al<^ 
constructed an SU(4) Majumdar-Ghosh model in a two-leg spin-3/2 ladder whose 
ground state is solvable exhibiting this plaquette state. Similarly to the quartetting 
order in the superfluid state, this plaquette order is of four sites without any site and 
bond spin orders. Upon doping, we suggest that a quartetting superfluid appears. 

Now let us speculate the phase diagram in the entire parameter regime as de- 
picted in Fig. [5] We expect that the long range Neel order can survive with a finite 
value of J2 as marked as phase C in Fig. [8l The plaquette ordered phase is gapped, 
thus it should be stabilized in an entire phase as marked phase C. It extends into a 
finite regime with J2 < Jo, and we also speculate that it survives in the entire region 
of Jo < J2. Furthermore, between phase C with the site spin order and phase A 
with the plaquette spin order, we suggest the existence of the magnetically ordered 
dimer phase B. The dimer made of two spin-3/2 particles possesses the internal 
spin structures. In phase B where J2 < Jo, the energy of the dimer singlet state 
is lower than the energy of the dimer quintet. Their superposition might give rise 
to an antiferromagnetically ordered spin-nematic dimer state. More analytic and 
numeric works are desired to confirm these speculations. 
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Fig. 9. ID phase diagram at half-filling (two particles per site) with t/rj,2 > (t is rescaled to 1). 
The phase boundary is marked with the solid line, and the 517(4) line (Uo = U2) is marked with 
the dashed line. A) The site singlet phase; B) The spin-Peierls (dimer) phase. 



6.4. spin-3/2 magnetism at half- filling 

At half- filling with Uo,U2 ^> t, the exchange model has been derived as Eq. fT9l 
In ID, the phase diagram of Eq. [19] can be obtained by applying the bosonization 
method directly to the Hubbard model after taking into account the 4fc/-Umklapp 
terms as 

H um ,4kf = Y r) R r)L + -^^rXl + h - c -> ( 55 ) 

where n R = ^ Ra R a /3^ Rl3 , XR = ^ Ra i Rva )ap^ R p, and tj l ,xl are defined corre- 
spondingly The bare values for A s>v are X s — Uoao,X v = L^ao where ao is the 
lattice constant. This term can be bosonized as 

HumAkf = - n/ 1 Tn cos(V47r0 c )((A„ + A s ) cos V^n^v + (An - A s ) cos Vin0 v 
2(ira) z I 

+ 2A.„(C0S \f^K(j)tl + COS %/47T</> t 2 ) j , (56) 

Together with Eq. [551 the RG equations at the half-filling can be calculated as 

rl\ 1 

(g c X v + g v X s + <±gtX v ). (57) 



dint 2?r 

Eq. [57|are solved numerically in the region Uo > 0, U2 > 0. Two stable fixed points 
are found as 

A : 9c = ~g v = 9t = -A s = X v -> +00, B : g c = g v = g t = X s = X v —> +00. (58) 
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with the corresponding phase diagram shown in Fig. [5J Phase A is the site singlet 
phase with the pinned values of the boson fields as (<fi c ) — (6 V ) — (<fiti) — {4>t2) = 0. 
All the charge orders and spin-Peierls orders vanish in this phase. Phase B is the 
dimerized phase with pinned boson fields as (<j> c ) — (<p v ) — {4m) — {<t>t2) — 0- As a 
result, the 2fc/ spin-Peierls (dimer) order is long range ordered. We notice that the 
SU(A) line of Uq = U2 is already inside the dimer phase, thus is not critical. This 
agrees with the previous RG result in Ref. The phase boundary is determined 
by numerically solving the above RG equations, which extends to the region with 
Uq < 1/2- The transition between phase A and B is Ising-like, which is driven by 
the competitions between cj> v and 9 V . 

The two dimension phase diagram of Eq. [19] is also challenging. As discussed in 
Sect. 13.31 at Uq < U2 Eq. [T9l reduces to the SO(5) rotor model. When the difference 
between U2 and Uq is small, or large compared to the inter-site exchange, the system 
is in the antiferromagnetic spin-nematic phase, or in the site singlet phase. In the 
regime of U2 < Uq, it is not clear whether the ground state still possesses magnetic 
long range order. An interesting numeric QMC result showed that along the 
SU (4) line, Eq. [TO] is in the gapless spin liquid phase, but more numeric work is 
desirable to confirm it. 

7. Summary 

In summary, we have reviewed the hidden 50(5) symmetry in spin-3/2 cold atomic 
systems. This symmetry is proved to be exact in the continuum model with s- 
wave interactions, and in the lattice Hubbard model in the optical lattices. This 
high symmetry provides a framework to understand various properties in spin- 
3/2 systems, including the Fermi liquid theory, mean field phase diagram, and 
the sign problem in QMC simulations, spin-3/2 systems support novel superfluid 
states including the quintet Cooper pairing state and the quartetting state. The 
topological defect of half-quantum vortex and the 5*0(4) Cheshire charge effect 
were discussed in the quintet superfluid. The existence of the quartetting state in 
ID systems and its competition with the Cooper pairing state were investigated. 
At last, we reviewed the magnetism in spin-3/2 systems showing many different 
features from the spin-^ systems. 

Taking into account the rapid progress in the cold atomic physics, we are opti- 
mistic that the spin-3/2 high spin cold atomic systems can be realized in experiment 
in the near future. We hope that the research on the hidden symmetry aspect can 
stimulate general interest in various properties in such systems. 
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